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Abstract:  The  standard  practice  of  burying  water 
and  sewer  lines  beneath  the  frost  line  in  cold  regions 
can  be  expensive  when  ledge  or  other  difficult  mate¬ 
rial  is  within  the  burial  depth.  If  the  pipeline  can  be 
buried  at  a  shallower  depth  and  still  be  protected 
from  freezing,  a  significant  savings  in  excavation 
costs  can  be  realized.  A  finite  element  (FE)  program 
was  developed  to  predict  frost  penetration  depth 
around  buried  utility  pipeiines.  The  program  was 


used  to  design  and  assess  the  feasibility  of  various 
insulation  configurations  around  a  water  line  buried 
within  the  frost-susceptible  depth  in  Berlin,  New 
Hampshire.  Extensive  temperature  monitoring  was 
performed  to  evaluate  both  the  insulation  design 
and  the  prediction  accuracy  of  the  FE  program.  The 
first-year  results  are  very  promising,  showing  good 
accuracy  between  the  FE  results  and  actual  tem¬ 
peratures. 
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Regions  Research  and  Engineering  Laboratory  (CRREL),  Hanover,  New  Hamp¬ 
shire. 
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vancement  Research  (CPAR)  program.  Tlie  report  was  technically  reviewed  by  Dr. 
Gary  E.  Phetteplace  of  CRREL  and  Dr.  John  Budinscak  of  U.C.  Industries,  Inc. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional 
purposes.  Citation  of  brand  names  does  not  constitute  an  official  endorsement  or 
approval  of  the  use  of  such  commercial  products. 
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Frost  Shielding  Protection  of  a 
Water  Line,  Berlin,  New  Hampshire 

BARRY  A.  COUTERMARSH 


INTRODUCTION 

In  areas  that  experience  freezing  temperatures 
for  any  appreciable  time,  protecting  water  lines 
from  freezing  is  a  major  consideration  in  the  de¬ 
sign  and  construction  of  water  systems.  The  nor¬ 
mal  procedure  is  to  bury  the  systems  below  the 
anticipated  frost  penetration  depth.  The  presence 
of  ledge  above  this  depth  can  greatly  increase  the 
cost  of  burying  the  system.  This  is  not  only  be¬ 
cause  of  the  expense  of  blasting  and  removing  the 
rock  but  also  because  frost  penetration  is  deeper 
in  areas  of  ledge,  requiring  an  even  deeper  pipe 
burial  depth  to  protect  against  freezing. 

Frost  shielding  is  the  practice  of  protecting  the 
water  line  from  freezing  by  adding  extruded  poly¬ 
styrene  insulation  to  retard  the  heat  loss  from  the 
water  within  the  line.  It  has  the  potential  of  re¬ 
ducing  the  burial  depth  of  the  system  from  what 
would  normally  be  necessary.  Balanced  against 
the  increased  cost  of  the  extruded  insulation  is 
the  decreased  depth  required  for  safe  burial  of 
the  pipe  and  thus  savings  in  both  the  volume  and 
cost  of  removed  material. 

The  major  impediments  to  routinely  insulating 
lines  are  the  lack  of  performance  data  and  design 
guidance  for  the  shield  configuration  and  insula¬ 
tion  thickness.  A  finite  element  (FE)  computer 
program  at  CRREL  can  model  two-dimensional 
frost  penetration  into  the  ground  and  describe 
both  numerically  and  visually  the  temperature 
regime  expected  in  and  around  any  potential 
shield  design.  This  FE  program  allows  the  de¬ 
signer  to  model  several  different  insulation  con¬ 
figurations  and  perform  "what-if"  types  of  calcu¬ 
lations  against  expected  temperature  conditions. 

This  report  illustrates  its  use  in  the  design  of 
an  insulation  shield  for  an  8-in.  (20-cm)-diam. 
water  line  in  Berlin,  New  Hampshire,  in  an  area 
where  ledge  is  present  to  the  surface.  It  details  the 


construction  of  the  shield,  methods  used  during 
construction  to  handle  the  2-  x  8-ft  (0.6-  x  2.4-m) 
extruded  polystyrene  insulation  boards  that  make 
up  the  shield,  and  the  temperature  sensor  layout 
we  used  to  monitor  the  shield's  performance.  The 
project  will  be  used  to  assess  the  performance  of 
both  the  FE  program  and  the  insulation  shield 
under  in-situ  conditions. 

Background 

This  frost-shielding  project  was  developed  un¬ 
der  the  Corps  of  Engineers,  Civil  Works,  Con¬ 
struction  Productivity  Advancement  Research 
(CPAR)  program.  Under  the  program,  the  City  of 
Berlin,  New  Hampshire,  Water  Works  and  U.C. 
Industries,  Inc.,  a  manufacturer  of  extruded  poly¬ 
styrene  insulation,  are  partners  with  CRREL  in 
demonstrating  the  concept  of  insulating  water 
lines  to  protect  them  from  freezing.  The  project 
has  essentially  two  components.  The  first  is  to 
develop  a  shield  design  using  a  CRREL-devel- 
oped  FE  program.  The  second  is  to  construct  and 
assess  the  performance  of  an  in-situ  shield  based 
upon  the  numerical  design.  More  information 
about  the  finite  element  program  is  available  in 
Coutermarsh  and  Phetteplace  (1991a,  1991b). 

The  practice  of  frost-shielding  utility  lines  has 
been  used  most  notably  in  Norway,  where  Per 
Gunderson  has  been  instrumental  in  its  develop¬ 
ment.  He  has  developed  nomographs  to  estimate 
the  shield  configuration  and  thickness  using  cer¬ 
tain  environmental  parameters.  The  Norwegians 
have  also  pioneered  collocating  sewer  lines  within 
the  same  shield  as  water  lines.  The  warmer  sewer 
lines  supply  heat  to  the  shield,  which  helps  to 
keep  inside-shield  temperatures  above  freezing 
(Gunderson  1975, 1989). 

It  was  the  success  of  Gunderson's  work  that 
encouraged  us  to  further  investigate  the  concept 
of  frost  shielding  in  this  country.  We  felt,  how- 


ever,  that  the  nomographs  were  of  limited  use  in 
conditions  where  several  different  soil  materials 
were  present,  the  geometry  of  the  pipeline  was 
nonstandard,  or  the  ambient  temperature  profiles 
were  erratic.  With  the  advent  of  more  powerful 
personal  computers,  we  felt  that  an  FE  program 
would  offer  the  designer  more  flexibility  in  choos¬ 
ing  shield  design  based  upon  the  conditions  pres¬ 
ent. 

The  city  of  Berlin,  New  Hampshire,  situated  in 
the  White  Mountain  range  of  northern  New 
Hampshire,  has  a  population  of  about  12,000.  The 
town's  geography  is  mountainous  with  large 
areas  of  shallow  ledge  throughout.  The  Berlin 
water  department  approached  CRREL  to  see  if 
we  could  offer  any  advice  on  the  freezing  prob¬ 
lems  they  were  experiencing  with  their  water  dis¬ 
tribution  system. 

The  water  in  the  distribution  system  is  from 
two  sources:  the  Ammonoosuc  Filter  Plant  and 
the  Androscoggin  Treatment  Plant.  The  Ammo¬ 
noosuc  plant  gets  its  water  off  the  bottom  of  a 
reservoir,  and  the  Androscoggin  plant  is  supplied 
from  the  Androscoggin  River.  The  Ammonoosuc 
plant  only  filters  and  chlorinates  the  water,  but 
the  Androscoggin  plant  provides  a  full  range  of 
treatment.  Ehiring  the  winter,  the  water  leaving 
the  Ammonoosuc  plant  tends  to  be  roughly  2°C 
warmer  than  the  water  leaving  the  Androscoggin 
plant. 

The  distribution  system  itself  is  old,  with  many 
of  the  lines  made  up  of  small  galvanized  2-in.  to 
6-in.  (5.1  to  15.2-cm)-diam.  pipes.  A  large  portion 
of  these  pipes  are  buried  in  or  on  top  of  ledge, 
tend  to  be  very  shallow,  and  are  prone  to  freez¬ 
ing.  The  frost  penetration  in  this  area  of  New 
Hampshire  can  be  1.5  m  (5  ft)  or  more.  The  prob¬ 
lem  is  so  acute  that  from  January  through  March 
the  water  users  on  lines  prone  to  freezing  are 
asked  to  run  their  water  continuously.  They  are 
given  an  abatement  on  their  water  fees  by  being 
billed  on  the  average  of  the  previous  three  months' 
usage  rather  than  on  their  actual  usage  during 
the  freezing  period.  This  causes  a  situation  where, 
during  the  winter  months,  water  is  being  run 
through  the  system  just  to  prevent  freezing.  The 
cost  of  this  water  is  in  effect  subsidized  by  all  the 
users.  The  excess  water  is  discarded  into  the  sewer 
system,  where  it  puts  a  greater  load  on  the  sewer 
treatment  system  and  titierefore  increases  treat¬ 
ment  cost. 

The  city  is  in  the  process  of  upgrading  its  wa¬ 
ter  system  to  correct  this  situation.  The  presence 
of  shallow  ledge,  however,  makes  this  an  ex¬ 


tremely  costly  procedure.  Ledge  excavation  is  rela¬ 
tively  expensive  because  of  the  need  to  blast  to 
break  up  the  material  for  removal.  The  situation 
is  exacerbated  because  frost  penetrates  much 
deeper  in  ledge  than  it  does  in  granular  soils.  This 
results  from  two  factors:  first,  the  ledge  has  a 
higher  thermal  conductivity  than  soil,  and  sec¬ 
ond,  because  it  has  no  moisture  content,  it  has  no 
latent  heat  (the  energy  used  when  changing  the 
phase  of  a  material,  e.g.,  changing  water  to  ice).  It 
can  be  seen  that  an  effective  shielding  method 
that  would  allow  a  shallower  burial  depth  than 
would  otherwise  be  necessary  would  have  the 
potential  to  save  a  substantial  amount  of  money 
and  resources.  U.C.  Industries,  Inc.  expressed  an 
interest  in  the  concept  of  frost-shielding  and  be¬ 
came  a  partner  with  the  water  department  to  pro¬ 
vide  guidance  and  insulation  for  the  project  as 
well  as  to  be  the  technology  transfer  point  of  con¬ 
tact  to  private  industry. 

This  report  describes  the  numerical  design  and 
construction  of  a  20.3-cm  (8-in.)  shielded  water 
line  in  Berlin,  New  Hampshire.  The  ongoing 
project  will  evaluate  its  performance  and  adjust 
the  design  procedure  accordingly  as  data  are 
evaluated. 


NUMERICAL  MODEL 


The  physical  dimensions  of  the  pipeline  and 
surrounding  material  are  used  to  construct  an  FE 
mesh  (Fig.  1).  Only  half  of  the  physical  configura¬ 
tion  needs  to  be  modeled,  since  the  other  half  will 
be  identical  to  the  modeled  half.  A  mesh  is  there¬ 
fore  constructed  that  represents  one  side  of  a  ver- 


Figure  1.  Example  of  a  finite  element  mesh. 
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tical  line  of  symmetry  through  the  insulation  and 
pipe.  This  mesh  consists  of  a  series  of  triangular 
elements  arranged  in  zones  that  accurately  repre¬ 
sent  the  physical  configuration  of  the  material 
being  studied.  The  density  of  the  element  spacing 
depends  upon  the  temperature  gradients  expected 
in  any  zone.  Areas  of  steep  gradients  (high  tem¬ 
perature  change)  should  have  elements  that  are 
spaced  closely  together.  Thus,  in  Figure  1  there 
are  areas  that  have  very  closely  spaced  elements 
representing  the  insulation  boards.  These  are  the 
darker  vertical  and  horizontal  bands  evident  in 
the  figure.  The  insulation  boards  are  placed  in  an 
inverted  U  immediately  around  the  pipe,  which 
is  the  darker  half-round  object  at  the  lefthand 
side  of  the  mesh.  The  dark  bands  extend  beyond 
the  boundaries  of  the  insulation  boards  because  it 
is  easier  to  generate  them  as  a  continuous  band 
throughout  the  mesh  and  then  just  designate  the 
board  itself  in  the  material  file  by  their  coordi¬ 
nates.  For  example,  the  top  horizontal  board  ex¬ 
tends  from  the  lefthand  side  of  the  mesh  to  the 


month  were  used  as  our  surface  temperatures  in 
the  model. 

Actual  surface  temperature  is  usually  different 
from  air  temperature  and  can  be  determined  by 
the  use  of  so-called  "n-factors."  These  factors  re¬ 
late  the  effect  that  surface  material  has  upon  the 
surface  temperature  (Lunardini  1981).  By  multi¬ 
plying  the  air  temperature  by  the  appropriate 
material  n-factor,  the  surface  temperature  can  be 
determined.  A  good  illustration  of  the  surface  tem¬ 
perature  difference  is  the  hot  temperatures  found 
in  summer  on  a  black  asphalt  surface,  which  will 
normally  be  several  degrees  hotter  than  the  air 
temperature. 

Using  air  temperature  as  our  mesh  surface 
boimdary  temperature  gives  us,  on  average,  a 
somewhat  colder  surface  than  the  actual  would 
be.  This  was  desirable  from  a  conservative  point 
of  view  so  our  shield  design  would  be  somewhat 
overdesigned.  Figure  2  is  a  graph  of  the  air  tem¬ 
peratures  determined  from  our  records.  The  time 
scale  is  in  months  and  represents  a  full  year  start¬ 


righthand  edge  of  the  vertical  band.  The  vertical 
board  extends  from  the  bottom  of  the  horizontal 
band  to  just  below  the  bottom  of  the  pipe.  Each  of 
the  other  materials  present  in  the  physical  con¬ 
figuration  is  designated  by  the  appropriate  coor¬ 
dinates  of  the  zone  of  elements  in  the  mesh.  The 
program  uses  this  information,  along  with  the 
thermal  conductivity,  density,  specific  heat,  la¬ 
tent  heat,  and  phase-change  temperature  of  the 
components  to  determine  heat  transfer  through¬ 
out  the  region  under  study.  More  information  on 
finite  element  (FE)  modeling  can  be  found  in  text¬ 
books  on  the  subject  (Segerlind  1984). 

Once  the  physical  configuration  is  determined, 
all  the  boundary  conditions  that  will  affect  the 
mesh  temperatures  have  to  be  determined  for  a 
full  year's  cycle.  In  our  case,  the  surface 


ing  in  January. 

Water  temperature  was  more  difficult  to  deter¬ 
mine.  We  obtained  monthly  treatment  and  filter 
plant  outflow  water  temperatures  for  1992  and 
part  of  1993  and  used  these  records  to  estimate 
the  water  temperature  in  the  pipe.  There  can  be  a 
significant  error  in  these  temperatures  because 
the  water  temperature  will  change  as  it  travels 
throughout  the  distribution  system.  How  it 
changes  is  not  known.  For  instance,  m  the  winter, 
the  water  from  the  plants  might  cool  significantly 
as  it  travels  through  pipes  that  are  buried  in 
ground  that  is  below  freezing.  Conversely,  since 
a  large  portion  of  Berlin's  water  is  from  a  surface 
source,  it  may  in  fact  be  cooler  than  the  groimd 
for  a  large  part  of  the  season.  The  preliminary 


temperature,  water  temperature,  and 
geothermal  heat  flux  (heat  generated 
within  the  core  of  the  earth  that  slowly 
filters  up  to  the  surface)  were  the  bound¬ 
ary  conditions  that  were  needed  for  our 
problem.  Along  the  vertical  line  of  sym¬ 
metry,  the  boundary  condition  is  adia¬ 
batic  (zero  heat  flux). 

The  surface  temperatures  used  in  the 
model  were  determined  from  monthly 
extreme  air  temperature  records  for  the 
Berlin  area  from  1926  to  1992.  From  the 


records,  1972  appeared  to  be  the  cold¬ 
est  year,  and  the  average  of  the  coldest  Figure  2.  Average  air  temperatures  from  1972  used  to  represent 
and  warmest  temperatures  within  each  FE  surface  temperature. 
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Figure  3.  Berlin,  New  Hampshire,  filter  plant  water  temperatures. 


data,  which  we  will  show  later,  indicates  that 
early  in  this  year  the  water  at  the  filter  plant  was 
colder  than  the  ground  for  portions  of  the  year. 
Figure  3  shows  die  water  temperatures  we  used 
in  our  simulation. 

In  our  numerical  simulations  we  model  to  a 
depth  of  10  m  (32.8  ft),  which  can  reasonably  be 
assumed  to  represent  the  depth  of  zero  annual 
temperature  amplitude  (Phukan  1985).  This  is  the 
depth  at  which  the  effect  of  the  changing  surface 
temperature  is  not  noticeable.  At  this  bottom 
boundary  we  apply  an  average  value  of  the  geo¬ 
thermal  heat  flux  of  0.063  W/m^  (Lunardini  1981) 
moving  upward  into  the  mesh. 

SECOND  AVENUE 

Figure  4  is  an  example  of  the  basic  FE  configu¬ 
ration  used  to  model  what  we  initially  expected 
the  soil  and  pipe  configuration  to  look  like  at  our 
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Figure  4.  Physical  layout  ofFE-modeled  domain. 


original  test  site  on  Second  Avenue  in 
Berlin.  The  pipe  is  20.3-cm  (8-in.)  I.D. 
ductile  iron  set  1.5  m  (5  ft)  deep.  Mate¬ 
rial  6  is  the  pipe  and  water,  material  4 
represents  the  backfill  area  around  the 
pipe,  material  7  is  the  extruded  poly¬ 
styrene  insulation  boards,  and  materi¬ 
als  2  and  11  are  the  surrounding  earth 
and  ledge  material,  respectively.  The 
specific  dimensions  and  materials  asso¬ 
ciated  with  each  zone  might  change,  de¬ 
pending  upon  what  the  designer  wishes 
to  experiment  with.  For  instance,  in  our 
modeling  we  varied  the  pipe  depth,  in¬ 
sulation  board  thickness,  and  the  width  of  the  top 
insulation  board  to  see  the  effect  these  changes 
would  have  upon  the  temperature  profiles. 

We  thought  that  the  pipe  on  Second  Avenue 
could  be  installed  about  1.5  m  (5  ft)  deep,  but 
ledge  was  possible  anywhere  from  1.4  m  to  2  m 
(4.5  ft  to  6.5  ft).  In  Figure  4,  material  11  was  ledge 
and  material  2  was  a  "silty  sand"  commonly  found 
in  the  area.  The  backfill  material  4  is  a  sand.  As 
mentioned  above,  we  varied  the  extruded  poly¬ 
styrene  insulation  (material  7)  dimensions  during 
the  modeling  process  to  investigate  the  effective¬ 
ness  of  different  shield  thicknesses  and  place¬ 
ment.  Table  1  lists  the  materials  and  their  prop¬ 
erties. 

During  the  numerical  simulations,  we  defined 
a  failure  as  the  0°C  isotherm  touching  the  pipe. 
This  is  a  conservative  approach,  since  water  will 
remain  at  0°C  for  some  time  before  it  freezes  due 
to  its  latent  heat.  For  the  initial  runs  we  did  not 
apply  any  water  temperature  boundary  condi¬ 
tion  to  the  pipe.  This  would  be  similar  to  a 
pipe  with  no  flowing  water.  Commonly,  wa¬ 
ter  flowing  in  the  pipe  brings  heat  into  the 
insulation  shield,  so  modeling  without  wa¬ 
ter  flow  is  usually  a  more  conservative  pro¬ 
cedure. 

The  output  from  a  numerical  simulation 
is  presented  in  graphical  form  by  either 
graphs  or  contour  plots.  By  watching  how 
the  0°C  isotherm  progresses  with  time,  the 
effectiveness  of  the  insulation  shield  can  be 
determined.  Figure  5  shows  a  time  series  pro¬ 
gression  of  five  isotherms  around  and  into 
a  shield.  Notice  how  the  isotherms  wrap 
around  the  shield  and  start  to  intrude  not 
only  at  the  top  but  also  at  the  lower  right 
comer.  This  is  due  to  the  higher-conductiv¬ 
ity  ledge  at  the  bottom  of  the  shield.  This 
contour  type  of  output  is  extremely  helpful 
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Table  1.  Numerical  model  material  data. 


Material  and  state 

Conductivity 

q/hrm 

Density 

(kglm^) 

Specific  heat 

(J/kg  x:) 

Latent  heat 
(]/n?) 

Material  4 

Sand 

Frozen 

Unfrozen 

5616 

6228 

1836 

1836 

784.8 

940.32 

45.3  X  106 

Material  6 

Water 

Frozen 

Unfrozen 

7984.8 

2170.8 

917.0 

998.2 

2096.2 

4183.92 

334  X 106 

Material  7 

Extruded  polystyrene 

104.04 

28.84 

1214.23 

Material  11 

Rock 

9242 

2700 

879 

Temperature  -2°C 
_  -1 


Figure  5.  FE  output  showing  isothermal  progression 
The  time  steps  are  12  h  apart. 
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in  assessing  the  shield  design.  A  large  amount  of 
information  can  be  displayed  in  a  format  that  is 
easily  interpreted. 

Our  previous  modeling  runs  and  the  earlier 
work  done  by  Gunderson  suggest  that  a  wide 
horizontal  layer  of  insulation  is  effective  at  slow¬ 
ing  frost  penetration.  We  modeled  this  on  the 
Second  Avenue  pipe  configuration  by  using  a  6- 
cm  (2.4-in.)-thick  layer  of  insulation  in  a  total 
width  of  2  m  (6.5  ft).  The  numerical  simulations 
predicted  that  this  would  prevent  the  0°C  iso¬ 
therm  from  contacting  the  pipe.  The  problem, 
however,  is  that  this  width  of  insulation  is  usu¬ 
ally  impractical  to  excavate  for  in  a  normal  utility 
line  installation.  That  is  why  it  is  often  preferable 
to  use  the  inverted  U  configuration  shown  in  Fig¬ 
ure  5. 

After  modeling  several  different  combinations 
of  insulation  thickness  and  shield  width,  using 
the  inverted  U,  we  found  that  a  10.2-cm  (4-in.)- 
thick  shield  with  a  total  minimum  width  of  101.6 
cm  (40  in.)  would  prevent  the  freezing  front  from 
touching  the  pipe. 

We  then  modeled  the  pipe  as  being  1.37  m  (4.5 
ft)  deep  and  found  that  there  was  not  a  practical 
insulation  thickness  that  could  prevent  the  freez¬ 
ing  front  from  touching  the  pipe.  The  next  ap¬ 
proach  we  took  was  to  see  if  water  flowing  in  the 
pipe  would  supply  enough  heat  to  prevent  freez¬ 
ing.  We  used  a  15.2-cm  (6-in.)-thick  shield,  1.22-m 
(4-ft)-wide,  and  ran  the  model  with  a  pipe  boimd- 
ary  temperature  to  simulate  water  flowing  in  it. 
We  then  stopped  the  flow  to  see  if  the  0°C  iso¬ 
therm  would  progress  to  the  pipe.  In  our  model 
this  is  a  two-step  procedure. 

The  first  step  of  the  numerical  simulation  is 
performed  with  a  specified  pipe  boimdary  tem¬ 
perature  derived  from  the  known  water  tempera¬ 
tures.  This  simulates  the  actual  temperatures  at 
any  time  we  might  need.  These  are  then  used  as 
starting  temperatures  in  the  second  step,  for  a 
numerical  simulation  without  the  specified  pipe 
boundary  temperatures.  As  noted  earlier,  this  is 
equivalent  to  a  situation  where  no  water  is  flow¬ 
ing  in  the  pipe.  We  then  run  the  model  for  any 
length  of  time  needed  to  determine  if  the  0°C 
isotherm  progresses  to  the  pipe. 

The  protection  of  the  water  pipe  in  the  above 
situation  is  dependent  upon  the  initial  tempera¬ 
ture  of  the  water,  the  time  of  year  the  water  stops 
flowing,  and  how  long  it  remains  off.  In  our  simu¬ 
lations,  we  showed  the  time  of  maximum  frost 
penetration  to  be  approximately  85  days  after  1 


January,  or  26  March,  at  about  2.5  m  (8.2  ft)  deep. 
If  the  water  flow  stopped  at  this  time,  our  model 
showed  it  to  be  safe  from  freezing.  We  then 
stopped  the  water  at  the  end  of  February  and  ran 
the  model  for  70  days  or  into  the  first  of  May.  In 
this  scenario  the  pipe  got  very  cold  but  did  re¬ 
main  (barely)  above  freezing. 

From  the  above  numerical  simulations,  we  felt 
that  a  1.2-m  (4-ft)-wide  shield  with  either  10.1-  or 
15.2-cm  (4-  or  6-in.)-thick  insulation,  depending 
upon  how  deep  the  pipe  was  buried,  would  be 
adequate. 

LABOSSIERE  STREET 
Initial  design 

When  the  initial  excavation  began  on  Second 
Avenue  it  was  discovered  that  the  existing  pipe 
leading  into  the  avenue  was  tightly  surrounded 
by  ledge.  It  would  have  been  very  difficult  and 
expensive  to  excavate  around  it  while  still  main¬ 
taining  water  service  to  the  dwellings  on  the  line. 
We  decided  to  move  to  our  second-choice  street — 
Labossiere.  The  water  line  in  this  street  was  a 
dead-end  line  about  123  m  (405  ft)  long,  and 
records  showed  that  the  existing  pipe  was  buried 
from  1.2  to  1.4  m  (4  to  4.5  ft)  deep.  We  dug  three 
test  pits  along  the  length  of  the  pipe,  and  these 
pits  indicated  that  the  closest  ledge  was  about 
1.37  to  1.52  m  (4.5  to  5  ft)  below  the  surface.  Since 
this  was  essentially  the  same  situation  we  had 
thought  we  had  at  Second  Avenue,  we  felt  that 
the  design  could  remain  the  same. 

Final  design 

It  was  quickly  discovered  during  the  initial 
excavation  for  the  new  pipeline  that  ledge  was  in 
fact  present  up  to  the  surface  for  nearly  the  entire 
length  of  the  pipeline.  In  an  incredible  stroke  of 
bad  luck,  the  only  three  places  that  showed  ledge 
down  to  1.4  to  1.5  m  (4.5  to  5  ft)  were  the  three 
spots  where  we  had  dug  our  test  pits.  Conse¬ 
quently,  we  found  ourselves  in  a  situation  where 
the  pipeline  was  being  installed  and  our  design 
was  suspect.  Again,  this  is  because  ledge  has  a 
much  higher  thermal  conductivity  and  contains 
less  moisture  than  soil.  These  two  conditions  al¬ 
low  the  freezing  front  to  advance  faster  and  deeper 
through  ledge  than  through  relatively  moist  soil. 

A  series  of  numerical  simulations  were  made 
that  modeled  ledge  to  the  surface  with  a  pipeline 
buried  1.5  m  (5  ft)  deep.  This  depth  was  close  to 
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where  the  pipeline  initially  started  out  at  the  be¬ 
ginning  of  the  street.  During  the  blasting  and  ex¬ 
cavation  we  tried  to  maintain  at  least  this  depth 
to  give  us  an  extra  buffer  against  the  cold  surface 
temperatures. 

The  numerical  model  predicted  that  the  pipe 
would  get  very  close  to  0°C  but  would  not  go 
below  it  by  using  our  final  design  from  Second 
Avenue:  a  1.5-cm  (6-in.)-thick  extruded  polysty¬ 
rene  insulation  with  a  total  shield  width  of  1.2  m 
(4  ft)  with  the  pipe  resting  at  the  1.5-m  (5-ft)-deep 
level. 

Several  points  should  be  made  about  our  de¬ 
signs  based  upon  the  numerical  simulations: 

1.  The  temperatures  used  at  the  model  surface 
are  air  temperatures  from  a  very  cold  year.  The 
temperatures  used  at  the  model  surface  should 
be  ground-surface  temperatures,  which  tend  to 
be  warmer  overall  than  the  air  temperatures  be¬ 
cause  of  the  effect  of  solar  energy.  This  study  will 
help  to  quantify  this  effect. 

2.  The  ledge  is  modeled  as  if  it  were  solid  un¬ 
broken  rock.  Natural  ledge  is  frequently  broken 
and  can  have  areas  of  water  flow  and  other  anoma¬ 
lies.  We  expect  the  effective  thermal  conductivity 
of  the  natural  ledge  to  be  less  than  our  modeled 
ledge. 

3.  Our  failure  criteria  could  be  much  too  con¬ 
servative.  We  "turn  off"  the  water  flow  and  as¬ 
sume  that  the  pipe  fails  if  afterwards  the  0°C  iso¬ 
therm  touches  it.  In  reality  there  is  probably  some 
water  flowing  frequently  within  the  pipe,  which 
would  bring  a  little  heat  into  the  shield  area. 

4.  A  great  uncertainty  exists  in  the  water  tem¬ 
perature  boundary  that  we  applied  to  the  pipe. 
As  mentioned  above,  we  used  the  water  tempera¬ 
tures  at  the  treatment  and  filter  plants  for  our 
pipe  boundary  temperatures.  This  may  or  may 
not  be  correct,  depending  upon  how  the  water 
temperature  changes  as  it  travels  through  the  dis¬ 
tribution  system. 

5.  Another  concern  in  this  study  is  the  effect 
that  Berlin's  very  cold  water  temperatures  will 
have  within  an  insulated  system.  Will  the  cold 
temperatures  quickly  remove  the  stored  heat 
within  the  insulation  shield,  much  faster  than 
would  have  occurred  had  the  heat  loss  only 
been  caused  by  heat  conduction  to  the  sur¬ 
rounding  ground?  One  of  the  ways  an  insu¬ 
lated  system  is  effective  is  by  preventing  the 
rapid  loss  of  this  stored  energy  to  the  sur¬ 
rounding  ground.  The  cold  pipe  tempera¬ 
tures  could  cancel  out  this  benefit. 


CONSTRUCTION 

General  physical  layout 

Labossiere  Street  is  a  dead-end  hillside  street 
on  the  northwest  side  of  Berlin.  Our  test  section 
started  at  the  intersection  of  Sixth  Avenue  and 
continued  up  the  roughly  17%  slope  of  Labossiere 
to  its  end  at  the  base  of  a  mountain  of  bare  rock 
that  projects  upwards  from  the  end  of  the  street. 

The  water  line  serving  the  users  on  the  street  is 
a  dead-end  line  that  nms  about  128  m  (420  ft)  up 
the  west  side  of  the  street. 

Pipe  construction 

The  shield  design  called  for  a  15.2-cm  (6-in.)- 
thick  layer  of  insulation  in  an  inverted  U  around 
the  pipe.  The  sides  of  the  U  were  61  cm  (2  ft)  high, 
with  the  bottom  of  the  legs  even  with  or  slightly 
below  the  bottom  of  the  pipe.  The  total  width  of 
the  shield  was  1.2  m  (4  ft).  The  shield  was  con¬ 
structed  of  5.1-cm  (2-in.)-thick,  4-  x  8-ft  SSE  boards 
of  Foamular®  250  extruded  polystyrene.  This 
material  meets  the  specification  requirements 
of  ASTM  C578,  Type  IV.  The  boards  have  lines 
scored  at  40.6,  61,  and  81.3  cm  (16-in.,  24-in.  and 
32-in.)  spacing  across  their  width,  which  make 
them  easy  to  break  at  the  job  site  and  leaves  a 
clean  even  line  where  they  break.  Consequently, 
the  61-cm  (2-ft)-high  sides  of  the  shield  could 
easily  be  obtained  from  the  4-  x  8-ft.  (1.2-  x  2.4-m) 
boards.  Table  2  lists  insulation  data. 

The  new  8-in.  ductile  iron  pipe  construction 
started  with  a  T  connection  into  the  existing  wa¬ 
ter  line  at  Sixth  Avenue  and  proceeded  up 
Labossiere  Street.  We  did  not  start  insulating  the 
line  from  this  T  because  of  the  presence  of  a  storm 
drain  line  that  ran  down  Sixth  Avenue  and  crossed 
over  the  new  water  line  just  up  from  this  connec¬ 
tion,  as  shown  in  Figure  6.  We  were  concerned 
that  cold  air  would  be  flowing  through  the 
underdrain  and  would  cool  down  the  water  pipe 
where  they  cross,  so  we  installed  a  thermocouple 
to  monitor  the  temperature  at  the  bottom  of  the 
drain  pipe.  The  stations  shown  in  Figure  7  start 
where  the  storm  drain  crosses  the  water  line.  The 


Table  2.  Insulation  data. 

Thermal  conductivity  Density  Specific  heat 
_ (J/hrm  ^C) _ (kg/m^)  (]/kg  °C) 

Foamular®  250  93.46  ^  28.84  1339.8 
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Figure  7,  Instrumentation  layout  along 
the  pipeline. 


insulation  shield  begins  about  2.4  m  (8  ft)  past  the 
underdrain. 

The  top  of  the  new  water  line  at  the  storm 
drain  is  about  1.5  m  (5  ft)  below  the  surface.  A 
backhoe  was  used  initially  to  excavate  for  the 
new  line,  as  the  ledge  that  was  encountered  was 
relatively  soft.  We  quickly  ran  into  very  compe¬ 
tent  ledge,  however,  and  the  ditch  line  had  to  be 
blasted;  a  larger  excavator  was  brought  in  to  re¬ 
move  the  debris.  We  tried  to  clear  the  rock  down 
to  1.8  m  (6  ft)  deep  with  a  minimum  1.2-m  (4-ft)- 
wide  ditch.  The  width  varied  depending  upon 
the  amount  of  rock  that  was  loosened  in  the  blast- 
ing. 

After  the  loose  rock  was  removed,  a  layer  of 
sand  about  15.2  cm  (6  in.)  deep  was  added  to 
cushion  the  pipe  on  the  ledge  at  the  bottom  of  the 
ditch.  We  positioned  some  insulation  jigs  made 
from  2-  X  4-in.-diam.  lumber  into  this  sand  and 
along  the  centerline  of  the  pipe  layout  to  hold  the 
bottom  of  the  side  pieces  of  insulation.  These  jigs 
prevented  the  sides  from  springing  out  at  the 
bottom  during  the  backfilling  operation.  At  the 
top  we  used  similar  jigs  to  firmly  hold  the  top  of 
the  sides  from  moving  in  or  out  during  the  back¬ 
filling.  These  top  jigs  also  maintained  a  1.2-m  (4- 
ft)  width  at  the  top  of  the  sides.  Figure  8  shows 
both  the  bottom  and  top  jigs. 


A  section  of  pipe  was  laid  on  top  of  these  jigs 
and  centered  in  the  ditch,  after  which  the  sides  of 
the  shield  were  installed.  It  is  important  to  pro¬ 
vide  even  mating  surfaces  and  to  stagger  the  in¬ 
sulation  joints  to  prevent  a  direct  path  for  heat 
loss  out  of  the  shield.  Three  layers  of  the  5.1-cm 
(2-in.)-thick  extruded  polystyrene  insulation  were 
necessary  to  get  our  final  design  thickness,  so  the 
boards  were  staggered  as  shown  in  Figure  9.  Any 
slope  changes  were  made  gradually  to  avoid  a 
gross  vertical  misalignment  of  the  side  boards. 
The  top  jigs  were  added  after  the  sides  were  in¬ 
stalled,  and  together  with  the  bottom  jigs  they 
held  the  boards  firmly  in  place  during  the  backfill¬ 
ing.  The  sand  backfill  was  alternated  between  the 
inside  and  outside  of  the  shield  sides  to  prevent 
the  walls  from  caving  in  under  the  soil  pressure. 
It  was  brought  up  even  with  the  top  of  the  side 
boards,  the  top  jigs  were  removed,  and  a  vibra¬ 
tory  compactor  was  used  to  compact  the  sand.  At 
this  point  the  sand  inside  the  shield  was  evened 
off,  usually  by  dragging  the  back  of  one  of  the  jigs 
across  the  tops  of  the  side  insulation.  This  made  a 
smooth,  level  surface  for  the  top  insulation  to  rest 
on.  The  top  insulation  was  staggered  in  the  same 
manner  as  the  sides  and  laid  down  over  the  pre¬ 
pared  surface  above  the  shield  sides.  Care  should 
be  taken  to  line  up  the  edges  with  the  sides  and  to 
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Figure  8.  Jigs  made  from  2-  x4-in.  lumber  to 
hold  the  insulation  during  installation. 


Figure  9.  Staggered  insulation 
board  joints. 


Figure  10.  Insulated  service  tap  entering  the  shield. 


avoid  any  gaps  that  would  provide  a  path  for 
heat  loss. 

Figure  10  shows  a  service  tap  entering  the 
shield.  These  taps  were  insulated  with  closed-cell 
pipe  insulation,  and  the  hole  made  in  the  insula¬ 
tion  shield  was  foamed  with  a  closed-cell  spray 
insulation.  Figure  11  shows  a  residential  sewer 


line  crossing  through  the  insulation  shield  at  ap¬ 
proximately  70  m  (231  ft)  up  the  pipeline.  The 
sides  of  the  shield  were  notched  so  the  insulation 
would  fit  over  the  pipe,  and  then  insulation  was 
fitted  into  these  notches  under  the  water  pipe 
after  the  sides  were  installed.  Some  spray  foam 
was  added  to  close  up  any  remaining  gap  be- 


Figure  11.  Residential  sewer  line  crossing  through  the  shield. 


tween  the  sewer  pipe  and  the  insulation,  as  shown 
in  Figure  11. 

The  shield  was  constructed  in  the  manner  de¬ 
scribed  above  up  to  the  end  of  the  line  at  approxi¬ 
mately  127  m  (417  ft).  It  would  have  been  desir¬ 
able  to  run  the  insulation  shield  1.5  to  3  m  (5  to  10 
ft)  past  the  end  of  the  pipe,  but  ledge  prevented 
this.  The  installation  of  the  shield  was  not  the 
limiting  factor  in  the  speed  of  the  pipe  laying;  the 
ledge  blasting  and  removal  took  considerably 
longer  than  the  time  needed  to  install  the  insula¬ 
tion.  Backfilling  was  slowed  slightly  by  the  extra 
care  taken  around  the  insulation  sides,  but  this 
time  was  minimized  after  we  gained  confidence 
in  the  ability  of  the  jigs  to  keep  the  insulation 
from  moving  during  the  process. 

INSTRUMENTATION  LAYOUT 

The  general  experimental  plan  called  for  moni¬ 
toring  temperatures  at  several  points  along  the 
pipe  with  thermocouples,  detailed  below,  as  well 
as  an  extensive  thermocouple  layout  at  a  test  sec¬ 
tion  about  73  m  (240  ft)  up  the  pipeline.  In  this 
application  we  used  copper/constantan,  type  T 
thermocouples.  The  temperature  data  is  automati¬ 
cally  recorded  by  a  Fluke  2286  datalogger  that 


was  installed  in  the  basement  of  a  house  on 
Labossiere  Street.  This  house  is  about  110  m  (360 
ft)  up  the  street  and  about  37  m  (120  ft)  from  our 
main  test  section.  We  also  installed  water  flow 
meters  in  the  pipe  at  two  locations  but  were  dis¬ 
appointed  to  discover  that  both  failed  to  function 
properly.  Figure  7  is  a  profile  view  of  the  as-built 
pipeline  with  several  points  of  interest  marked. 
Three  immersion  thermocouples  were  installed 
into  the  pipe  at  3.1,  73,  and  125  m  (10,  240,  and 
410  ft)  from  the  beginning  of  the  insulation.  These 
allow  us  to  monitor  the  water  temperature  as  it 
progresses  up  the  street.  At  our  main  test  site  we 
have  also  installed  thermocouples  around  the  out¬ 
side  of  the  pipe  near  the  middle  immersion  ther¬ 
mocouple.  This  will  allow  us  to  tell  if  there  is  any 
significant  difference  between  the  temperatures 
of  the  water  and  the  outside  of  the  pipe. 

The  bottom  immersion  thermocouple  gives  us 
the  temperature  of  the  water  near  the  entrance  of 
the  test  pipe.  This  will  be  a  baseline  temperature 
of  the  water  as  it  enters  the  test  section,  before 
any  influence  of  the  shield.  A  flow  meter  was 
installed  at  6.1  m  (20  ft)  up  the  pipe  but,  as  men¬ 
tioned  before,  it  is  not  working  properly.  We  had 
hoped  to  know  the  water  flow  into  the  test  sec¬ 
tion  by  this  meter.  Another  thermocouple  was 
placed  on  the  outside  of  the  pipe  at  about  24.4  m 
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Figure  12.  Sewer  line  crossing  with  thermocouple 
layout. 


file  in  Figure  13.  We  also  installed  an  immersion 
thermocouple  and  a  flow  meter  at  this  location. 
As  before,  this  flow  meter  does  not  appear  to  be 
working.  The  thermocouple  layout  at  the  main 
test  section  will  enable  us  to  characterize  the  ther¬ 
mal  environment  around  the  shield.  The  control 
string  will  give  us  temperatures  that  we  would 
expect  to  find  in  an  undisturbed  area.  The  tem¬ 
peratures  obtained  at  the  surface  of  the  test  sec¬ 
tion  will  be  used  as  our  surface  boundary  tem¬ 
peratures  in  the  numerical  model,  and  the 
temperatures  from  the  immersion  thermocouple 
will  be  used  as  our  pipe  boundary  temperatures. 

Three  more  thermocouples  were  installed  at 
the  site.  One  is  on  the  normally  dry  portion  of  a 
hydrant  standpipe  just  above  the  main  test  sec¬ 
tion  at  about  104  m  (340  ft).  There  is  an  immer¬ 
sion  thermocouple  at  the  end  of  the  pipeline  at 
125  m  (410  ft)  and  an  outside  air  temperature 
thermocouple  off  the  outside  north  wall  at  73 
Labossiere  Street.  We  planned  to  compare  air  and 
surface  temperatures  using  this  thermocouple  and 
the  previously  mentioned  surface  thermocouples 
at  the  main  test  section. 


(80  ft).  The  sewer  line  crossing  at  70.4  m  (231  ft) 
gave  us  an  opportunity  to  measure  the  effect,  if 
any,  that  the  relatively  warmer  sewer  line  would 
have  on  the  temperatures  within  the  shield.  Fig¬ 
ure  12  shows  the  sewer  line  and  the  placement  of 
the  four  thermocouples  at  this  location. 

The  main  instrumentation  section  was  at  ap¬ 
proximately  73  m  (240  ft)  along  the  pipe.  Here  we 
installed  40  thermocouples,  as  shown  in  the  pro- 


PRELIMINARY  DATA 

Air  temperatures 
and  frost  penetration 
The  winter  of  1994-95  was  unusu¬ 
ally  warm  and  was  therefore  not  an 
optimum  test  for  the  shielding  con¬ 
cept.  Some  of  the  temperatures  col¬ 
lected  from  the  beginning  of  data  col¬ 
lection  on  9  September  1994  imtil  17 
March  1995  are  shown  in  Figure  14, 
Here  the  outside  air  temperature  is 
compared  with  the  numerical  sur¬ 
face  boundary  temperature  used  in 
the  design  model.  As  mentioned  pre¬ 
viously,  this  FE  boundary  tempera¬ 
ture  was  determined  from  air  tem¬ 
perature  records  in  the  Berlin  area.  It 
can  be  seen  that  the  numerical  tem¬ 
peratures  were  generally  colder  than  the  air  tem¬ 
peratures.  Probably  just  as  important  is  that  when 
we  did  experience  cold  temperatures,  they  did 
not  last  very  long.  Notice  that  the  temperatures 
above  0°C  tend  to  last  longer  than  those  below 
O'^C,  and  there  is  never  a  prolonged  cold  spell. 
Under  this  type  of  temperature  profile  the  ground 
never  has  a  chance  to  maintain  any  appreciable 
frost  penetration.  In  our  situation  in  ledge,  the 
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Figure  14.  Air  temperature  vs.  FE  surface  temperature. 


frost  penetration  reacts  quickly  to  surface 
temperature  conditions.  Figure  15  shows 
how  rapidly  the  0°C  isotherm  changed 
depth  in  the  ledge,  depending  upon  the  sur¬ 
face  temperatures.  The  maximum  frost  pen¬ 
etration  that  we  recorded  on  Labossiere 
Street  was  about  1.1  m  (3.5  ft). 

Pipe  temperatures 

Berlin  has  very  cold  water  in  its  distribu¬ 
tion  system,  which  makes  it  a  challenge  to 
design  effective  insulation.  Figure  16  is  a 
graph  of  water  temperatures  at  the  treat¬ 
ment  plant  and  at  the  immersion  thermo¬ 
couples  in  the  pipe  at  Labossiere  Street.  It  is 
evident  from  this  figure  that  in  September 
the  water  at  the  plant  is  warmer  than  at  the 
beginning  of  Labossiere  Street  (the  lower 
immersion  thermocouple)  and  has 
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Figure  15.  Control  string  vertical  temperature  profiles. 


therefore  lost  heat  to  the  groimd  as  it 
flowed  to  Labossiere.  As  the  winter 
season  progresses,  there  are  periods 
where  the  treatment  plant  water  has 
warmed  up  as  it  has  traveled  to 
Labossiere. 

The  lower  and  middle  immersion 
thermocouples  are  roughly  the  same 
depth  below  the  surface,  and  the  up¬ 
per  thermocouple  is  1.2  to  1.5  m  (4  to 
5  ft)  deeper.  We  can  compare  the  low¬ 
er  and  middle  thermocouples  (Fig. 
17)  and  see  that  the  middle  thermo¬ 
couple  is  about  0.5°C  warmer  than 
the  lower,  showing  the  temperature 


rise  of  the  water  as  it  travels  up  the 


street.  The  upper  thermocouple  shows  the 


warmest  temperatures,  but  this  would  be  ex¬ 
pected  because  of  its  greater  depth.  This  makes 
it  difficult  to  use  the  upper  thermocouple  to 
draw  any  conclusions  about  any  beneficial 
effect  from  the  insulation. 

We  were  concerned  that  cold  air  flowing 
through  the  storm  drain  that  passes  over  the 
water  pipe  at  the  beginning  of  Labossiere  would 
hasten  freezing  in  the  immediate  area.  As  can 
be  seen  from  Figure  18,  the  temperatures 
recorded  at  the  bottom  of  the  storm  drain  were 
in  fact  warmer  than  the  lower  immersion  ther¬ 
mocouple  nearby.  We  do  not  know  if  this  would 
hold  in  a  year  with  colder  air  temperatures. 

Shield  temperatures 

Figure  16.  Comparison  of  lower  and  middle  immersion  In  comparing  temperatures  inside  and  out- 


thermocouples  and  filter  plant  temperatures. 


side  the  shield  at  the  1.5-m  (5-ft)  depth  (Fig. 


Figure  17.  Immersion  temperature  vs.  ground  temperature.  Figure  19.  Inside  vs.  outside  shield  temperatures. 
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Figure  18.  Comparison  of  lower  immersion  thermocouple,  Figure  20.  Pipe  vs.  ground  temperature, 

underdrain,  and  earth  temperatures. 


19),  it  is  evident  that  the  inside  temperature  is  NUMERICAL  MODEL 
colder  than  the  adjacent  outside  temperature.  This  ASSESSMENT 

is  a  result  of  the  cold  water  temperatures  in  the 

shield  and  the  fact  that  it  was  a  warm  year.  In  a  The  first  shield  was  designed  using  air  tern- 

normal  year  we  would  expect  the  inside  shield  peratures  that  we  hoped  would  be  representative 

temperature  eventually  to  become  warmer  than  of  an  extremely  cold  year.  This  was  a  conserva- 

the  outside  when  the  colder  temperatures  pen-  tive  approach  both  from  the  climatological  stand- 

etrate  deeper  into  the  groimd  as  the  year  pro-  point  as  well  as  the  numerical  modeling  stand- 

gresses.  Figure  19  also  shows  the  control  tern-  point  because  the  numerical  model  uses  surface 

peratures,  which  are  in  ledge,  becoming  colder  temperature,  which  is  usually  warmer  than  air 

than  the  equivalent  depth  at  the  pipe,  which  is  in  temperature.  Figure  21  shows  the  actual  surface 

sand  backfill.  It  is  evident  in  Figure  20  that  the  temperature  from  our  control  string  on  Labossiere 

ground  temperature  at  the  depth  of  the  shield  did  and  two  regression  lines.  One  regression  repre- 

not  get  colder  than  the  water  temperatures.  sents  the  FE  boundary  temperature  used  in  our 
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Figure  21.  FE  surface  temperature  vs.  actual  surface  tem¬ 
peratures,  Laboissiere  Street. 


peratures  from  the  same  time  period  were 
used  as  our  pipe  temperatures. 

Figures  19  and  20  also  show  the  results 
from  our  numerical  model.  In  Figure  19,  the 
numerical  model  overpredicts  the  inside 
shield  temperatures  from  September  until 
about  the  first  of  December,  but  it  does  very 
well  until  the  end  of  data  in  March.  The 
overprediction  is  in  part  because  our  nu¬ 
merical  pipe  temperatures  (from  the  treat¬ 
ment  plant)  are  warmer  than  the  actual  pipe 
temperatures  during  September  to  Decem¬ 
ber.  The  outside-the-shield  numerical  tem¬ 
peratures  are  at  first  warmer  than  the  ac¬ 
tual,  until  about  mid-November.  From  this 
point  until  the  end  of  December  they  repre¬ 
sent  the  actual  temperatures  well.  From  then 
until  the  end  of  data  in  March  our  numeri¬ 
cal  model  is  somewhat  colder  than  the  ac¬ 


tual  temperatures. 

The  numerical  model  did  poorly  on  pre¬ 
dicting  the  overall  frost  penetration  as  com¬ 
pared  with  the  control  string  temperatures 
in  ledge.  The  maximum  frost  penetration 
predicted  by  our  numerical  model  was 
roughly  1.8  m  (5.9  ft)  in  about  the  mid  to 
end  of  February.  The  actual  maximum  frost 
penetration  was  about  1.1  m  (3.5  ft)  during 
the  first  part  of  March. 

The  colder  numerical  temperatures  and 
deeper  modeled  frost  penetration  could  be 
due  to  our  ledge  thermal  conductivity  and/or 
the  nature  of  the  regression  model  and  how 
12  Nov  2  Dec^  22  Dec  1 1  Jan  31  20  Feb  12  Mar  ^  represents  the  actual  temperatures. 

I  Figure  22  shows  a  portion  of  the  control 

Figure  22.  FE  surface  vs.  actual  temperatures.  surface  temperatures  from  the  beginning  of 

November  until  the  beginning  of  March. 
Notice  that  once  the  regressed  surface  tern- 
shield  design,  and  the  other  is  one  that  fits  the  perature  line  goes  below  0°C  it  remains  there  un¬ 
actual  surface  temperature.  When  the  two  regres-  til  spring.  In  a  normally  cold  year  the  actual  tem- 

sions  are  compared  it  shows  that  the  actual  sur-  peratures  are  usually  modeled  well  by  this  type 

face  temperature  was,  on  average,  over  5.5°C  of  regression.  However,  in  the  winter  of  1994-95, 

warmer  than  the  numerical  surface  temperature.  the  actual  temperatures  rose  above  this  regressed 

To  better  determine  how  the  numerical  model  line  and  above  freezing  for  a  significant  amount 

performed,  the  top  and  pipe  boundary  conditions  of  time,  as  shown  in  the  figure, 

were  changed  to  reflect  the  actual  temperatures  To  check  the  effect  that  the  regression  line  has 

more  closely.  To  do  this,  we  need  a  year's  surface  upon  frost  penetration,  we  used  actual  surface 

and  pipe  temperatures  covering  the  model  year  temperatures  rather  than  the  regressed  tempera- 

of  interest.  We  had  been  gathering  data  on  Sec-  tures  as  the  top  boundary  condition.  The  0°C  iso- 

ond  Avenue  for  over  a  year,  and  these  surface  therm  reacted  much  more  realistically  in  this 

temperatures  were  compared  with  the  available  scenario,  with  the  ground  temperatures  rising 

Labossiere  temperatures  and  adjusted  to  repre-  above  0°C  during  some  of  the  warmer  surface 

sent  the  Labossiere  temperatures  from  March  1994  temperature  times.  The  maximum  frost  pene- 

to  March  1995.  The  treatment  plant  water  tern-  tration  time  was  still  around  the  mid  to  end  of 
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February  and  decreased  to  1.7  m  (5,6  ft),  a  small 
improvement. 

As  mentioned  earlier,  we  model  ledge  as  a  con¬ 
tinuous  mass  of  solid  rock  where  actual  ledge 
might  not  be  so.  During  the  construction  of  the 
pipeline  it  was  evident  that  portions  of  the  ledge 
were  fractured,  with  seams  of  water  within  some 
of  the  fractures.  These  areas  would  greatly  influ¬ 
ence  the  thermal  conductivity  of  the  area  and 
would  slow  frost  penetration  compared  to  solid 
rock.  The  above  analysis  suggests  that  we  have  a 
higher  numerical  ledge  thermal  conductivity 
than  is  appropriate. 

PRELIMINARY  CONCLUSIONS 
AND  RECOMMENDATIONS 

The  winter  of  1994-95  was  too  warm  to  pro¬ 
vide  us  with  a  good  test  of  the  shielding  concept. 
We  still  do  not  know  what  the  water  tempera¬ 
tures  will  be  at  the  test  site  during  a  very  cold 
year  or  what  effect  they  will  have  upon  the  shield's 
performance. 

It  appeared  that  the  shield  was  providing  some 
protection  to  the  water  line,  as  the  temperature  of 
the  water  increased  as  it  flowed  up  the  hill.  It  is 
difficult  to  draw  any  firm  conclusions,  however, 
since  the  ground  at  the  pipe  depth  remained 
above  freezing  throughout  the  winter. 

Another  area  of  concern  is  the  performance 
of  the  service  lines  to  the  residences.  These 
small-diameter  lines  are  relatively  shallow  and 
uninsulated  and  would  seem  to  be  susceptible 
to  freezing. 

The  numerical  model  appears  to  be  function¬ 
ing  well  where  we  have  accurate  input  to  check  it 


against.  In  our  future  comparisons  we  will  be 
using  the  actual  surface  temperature  as  our  top 
boundary  condition  to  see  if  we  can  more  closely 
model  the  maximum  frost  penetration  depth.  The 
results  of  these  runs  will  help  tell  us  if  the  ledge 
thermal  conductivity  needs  to  be  adjusted  and 
what  effect  inaccuracies  in  the  material  thermal 
conductivities  have  upon  our  model  results, 
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